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Abstract 



We investigate dynamics of overlapping vortices in the nonlinear Schrodinger equation, the nonlinear heat equation and in 
the equation with an intermediate Schrodinger-difFusion dynamics. Because of formal similarity on a perturbative level we discuss 
also the nonlinear wave equation (Goldstone model). Special solutions are found like vortex helices, double-helices and braids, 
breather states and vortex mouths. A pair of vortices in the Goldstone model scatters by the right angle in the head-on collision. 
It is found that in a dissipative system there is a characteristic lenght scale above which vortices can be entangled but below 
'which the entanglement is unstable. 

> ■ 

PO Introduction 



dynamics of the superfluid helium condensate at zero temperature U was proposed to be described by the nonhnear Schrodinger 
'eg^tion . Since then many approaches have been developed to study the dynamics of vortices being topological excitations of 
,tlj(e-jNLSE. Lund |Q derives equations of motion for 3-dimensional vortex filaments with the help of the effective action method. 

effective model takes the form of an action at a distance theory. As anticipated such a method can not describe many 
do^Slils of the dynamics. Lee Q rederives Lund's equations but in addition takes into account sound radiation from a moving vortex 
filMient. His method originates from earlier similar developements in relativistic field theories 

GThe dynamics of well separated vortices on a plane has been studied by Neu . The basic assumption of his method is that 
tb^modulus of the scalar field around a given vortex is undisturbed by other distant vortices. Interactions are mediated only by 
tl y]p hase of the complex scalar field. Such a simplification is not reliable for vortices with overlapping cores. 

Q[n this paper we are going to consider the dynamics of vortices when their cores strongly overlap. The method can be outlined 
as^Uows. We take as a background configuration a planar vortex solution with winding number say 2. Then we consider small 
fliiptuations around this background. The double vortex can be viewed as a superposition of two unit vortices. As the double vortex 
!^higher energy then two single vortices, it is likely to split into vortices with winding number 1. This splitting is described in 
nperturbative calculation as a mode with dipole symmetry. Its energy is negative as it should be for a splitting mode. One can 
algt^ consider axially symmetric modes which do not split the double vortex but instead they change its width. They are usually 
called breathers or pulsons. Another possibility is a mode with Z3 symmetry. The double vortex might uniformly split into 3 unit 
vortices leaving one antivortex at the center. The perturbative calculation shows this way of splitting does not take place for the 
double vortex. However for a vortex with winding number 3 there is a Z4-symmetric mode which describes a decay of the triple 
vortex into 4 unit vortices with one antivortex left at the origin. We generalize these considerations allowing the modes to vary 
along the line of the background vortex. It this setting the double vortex can be split into parallel unit vortices and later on these 
vortices can be twisted to form a twisted pair or a double helix. It turns out that there is a critical wave-lenght at which such a 
double-helix is static. 

Our analysis is carried out not only for the nonlinear Schrodinger equation (NLSE) but also for a diffusive model we call after 
Neu 1^ the nonlinear heat equation (NLHE). We discuss an intermediate Schrodinger-diffusive dynamics case and that of the 
nonlinear wave equation (NLWE) also known as the Goldstone model. The physical reason to consider diffusive and hamiltonian 
dynamics in one paper is that systems like superfluid helium change their dynamical properties according to external conditions 
such as temperature. At zero temperature the dynamics is hamiltonian but as we move towards the phase transition it gradually 
becomes more diffusive. Thus the different equations can describe one physical system but in different regimes. There are also 
important mathematical reasons for the unified treatement. First of all there is a set of static solutions, like static double helices, 
which are common to all the considered equations. The eigenvalue problem for time-dependent modes turns out to be formally the 
same for the NLWE and the NLHE. Solving this problem we kill two birds with one stone. The other reason is that there is a direct 
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corespondence between the modes in the NLHE case and those in the NLSE. Although the at first sight temptating analytical 
continution in the complex time plane fails, a large class of perturbative NLHE's solutions can be mapped into corresponding 
perturbative solutions of the NLSE. The analytical continuation does not apply to the modes themselves but it does apply to their 
eigenvalues. Some NLHE's modes do not map to the NLSE. However they still can be mapped to the intermediate diffusive equation 
with some admixture of the Schrodinger dynamics. 



2 Field equations and vortex solutions 

We will consider three types of nonlinear field equations. The first two are the nonlinear heat equation (NLHE) 



V',t = A7^+(f-|7^nv>^-^ (f) 



and the nonlinear Schrodinger equation (NLSE) 



-i7^,t = A^ + (l-| . (2) 

where A = is a Laplacian. The free energy and the Hamiltonian are 

F = H = jd^x [VrV^ + 1(1- I I)'] . (3) 
The third equation is the nonlinear wave equation (NLWE) 

V',tt-AV'-(l-|VnV'=^ = , (4) 

with the Lagrangian 

L = / d'x[d,rd^^-l{i-\^m . (5) 
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For time-independent fields all the three models reduce to the same static equation 

AV+(1- I V l^)V' = . (6) 

It is well known that such an equation admits topological vortices. A rotationally symmetric vortex solution can be obtained in the 
form 

^(r,0)=/„(r)e™« , (7) 
where {r,9) are polar coordinates. Substitution of this Ansatz to Eq.(^ yields 

A„/„ + (1 - /2)/„ = , (8) 

where A„ = + 7^ — 72"- The boundary condition at infinity is /„(oo) = I. The phase of the complex scalar field is multivalued 
at the origin so its modulus must vanish there, /n(0) = 0. If we restrict to regular solutions the asymptotes at the origin and at 
infinity will turn out to be 

Mr) « /or" + ... , r^O , 

TL 

From now on we will restrict to positive n. 



3 Fluctuations around vortex background in nonlinear heat and nonlinear wave 
equation 

Let the background solution be a straight-linear vortex along the z-axis with winding number n and the wavefunction -0, see (0..^. 
The particle density of the background will be denoted by p — The equation satisfied by the field fluctuation d^j} within the 
NLHE is 

- ^V-.t = -A,5V^ + (2p - 1)(5V + pe^"^5V* ■ (10) 
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One could try decomposition of the solution into Fourier modes both in 9 and z, say e*''*+*'^^. However such a single Fourier mode 
can not be a solution to Eq.([l0|) because of the last term, pe^"^^ Si/j* , which couples Fourier modes into pairs. The general solution 
can be obtained as a sum (or an integral) over pairs of Fourier modes like 

^E, fc,p)/oe-^*[e**'-^e*P%i(r) + e-*'=^e^(2"-P)%2 W] (11) 

or 

e(i?,fc,p)/oe-^*e-'=-[e*f«ui(r)+e*(2„-rt«^^(^)] ^ (^2) 

where f's are complex amplitudes of the modes and /o is a constant from Eqs.(^ which is introduced to provide the modes with 
a convenient normalisation. The z-dependence factorizes out because the background solution is z-independent. p's are integers 
which can range from —oo to n. Modes with p > n are identical to those with p < n e.g. the p = n+l modes can be identified with 
the p ^ n — 1 modes. The eigenvalues (E, k) can in principle take arbitrary real values but only the modes with finite energy (or 
free energy) are physical. The functions ui,M2 in Eq.(^) are solutions to the eigenvalue problem 

[E — k'^)ui = — ApMi + (2/9 — l)ui + pu2 , 

{E - k'^)u2 = -A2n-pU2 + i2p - 1)U2 + pui , (13) 



where A's are defined below Eq.(|8D. The functions ui,U2 in Eq.(|l2|) are solutions of the same eigenvalue problem as in Eqs.(|lJ 
but with k^ replaced by — fc^. From now on we will allow k^ to be a positive or negative real number, fc^ > will refer to the modes 
in Eq.(^T|). fc^ < will refer to those in Eq.([l^) so that k is to be replaced by Im{k) there. 
We can simplify Eqs.(|l3|) by a replacement E — k"^ = 

ujui = -ApUi + {2p - l)ui + pu2 , 

WU2 = -A2ri-pU2 + (2jO - 1)W2 + pUl . (14) 

With equations ( |l^ at hand we are able to give a physical interpretation to the modes with various p. If one restricts to regular 
solutions one can easily find the asymptotes of the profile functions close to the origin 

Ui{r) « M^r'^' , 

U2{r) w M°rl2«-Pl + ... , (15) 

where u\, U2 are real constants which have to be choosen so as to meet the conditions of regularity and fast convergence at infinity. 
Eqs.(p^ are a set of linear homogenous differential equations so the constants can be multiplied by a common factor. In other 
words we can redefine the constants by an appropriate rescaling of the overall amplitudes ^ in Eqs.(ll;l^). Thus, provided that 



u\ 7^ 0, we can choose u\ ~ 1 and call the still free = a. In this case we are left with just one coefficient a which we can vary 
to obtain solutions with acceptable asymptotes at infinity. For = the coefficient itj can be rescaled to 1, = 1. In this case, 
which is a set of measure zero, there is no free parameter to remove singularities or slowly convergent asymptotes at infinity and 
it turns out that there indeed is no solution. Thus from now on we restrict our attention to the first case 

Mi(r) « rl^l + ... , 

M2(r)«arl2"-?'l+... . (16) 
In the limit of very small w = x + iy or close to the vortex axis, for the modes in Eq.(|ll|), the total scalar field looks like 

V; = + 5V ~ li'" + ^{E, k,p)e-^'+''''wP (17) 

for p> 0. The zeros of this polynomial coincide with the zeros of the total scalar field. 

Let us consider the double vortex background, n = 2. For p = we have ip ^ w'^ + ^{E, k, 0)e~^*+*'^^. At t = and on the plane 
z = we have two zeros at complex roots of — ^. For E > these zeros will with time shrink down to w = 0. For nonzero k the two 
lines of vanishing scalar field form a double helix. If there is a solution for k then there is also a solution for —k. We can combine 
the two to obtain e.g. ip ~ + ^e~^* cosfcz with ^, say, real and positive constant. Let us fix time t but consider what happens 
as we vary z. Close to z = the zeros lie on y-axis. As we increase z the zeros shrink to w — but then reappear but rotated by 
the right angle. This kind of solution will be called a double-braid. 

For general n,p we get n — p uniformly split vortices and p vortices left at the origin. If p < there are —p antivortices at the 
origin. The case n = p is exceptional. There is no splitting but just a change of vortex width. 

Similar mode decomposition as for the NLHE (0,n2h can be performed for the NLWE, 



C(£;,fc,p)/o[e''^*+*'="ui(r)+e-*^*-*'="e*(2"-f)%2(r)] . (18) 
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Substitution of this Ansatz to the NLWE hnearised in fluctuations 

^2 

- ^ -^6ij+{2p-\)5iP^ pe^^'^^i^* . (19) 

once again leads to Eqs.(|l^) but this time lu = — k^. The Ansatz ( [l8| ) can be generalised to admit or negative. If say E^ 
is negative one should replace the two time-dependent exponents in the Ansatz (^8|) by one exponent e^-f™(^)* standing in front of 
the mode. An analogous rearrangement should be done if k'^ < 0. 

Since both the fluctuations in the NLHE and in the NLWE satisfy Eqs.(|lJ) it is the highest time to look for their solutions. 

3.1 Asymptotes of the splitting modes at infinity 

At infinity the asymptotic form of Eqs.(|l^) is 



These equations can be combined to give 



= Aui + {UJ ~ l)ui — 1*2 , 

= AU2 -Ul + {u}- 1)U2 . (20) 



Au+ + (w - 2)u+ , 

= Au_ + wit^ , (21) 



where m+ = ui + U2 and u_ = ui — U2. 

For UJ > 2 both u+ and u_ approach combinations of Bessel functions. In the range < a; < 2, u+'s asymptote is a linear 
combination of the modified Bessel functions while m_ still falls down like l/\/r. Finally for < both and u_ are independent 
combinations of the modified Bessel functions. The cases oi ui — 0,2 need a separate treatement. 

We can accept only localised modes with finite energy or free energy. With such a restriction we have to exclude both expo- 
nentially divergent solutions and all the solutions with Bessel function-like asymptotes. For ui = 2 there are slowly falling down 
oscillations in the asymptote of u_. Thus the range of w where we can look for acceptable solutions is restricted to w < 0. In this 
range there are two divergent asymptotes at infinity. For any given uj we can use the parameter a in Eq.(|l6|) to fine tune one of 
the divergent asymptotes to zero. The second can be removed by an appropriate choice of lo. Such a distinguished to is just the 
eigenvalue we are looking for. 

3.2 Bound states for p = 0..n — 1 

There are analytical solutions for p = n — 1 and w = 0. They can be constructed as 5ip = dx'il^{r,9), where ■0 = /(r)e*"* is the 
background n-vortex solution. The profile functions of such zero modes are 



2'-' ^ ' r ' 

u.^lifir)-!^] . (22) 
2 r 

The energy of these zero modes similarly as the energy of the background vortex is logarythmically divergent. They give rise to the 
following solution of the NLWE for small t 

5i^/fo^-vt[u,ir)+e^''u2{r)] , (23) 

where t; is a complex velocity. For n = 1 the zero of the scalar field coincides with the zero of the polynomial w — vt. The solution is 
just a perturbative approximation to a planar vortex moving with a constant velocity. This solution becomes more interesting if we 
admit k^ to excite travelling waves with a dispersion relation E"^ — k^ = 0. The line of vanishing scalar field coincides with zeros 
of a polynomial z) — Fi(t + z) ~ F2(t — z), where -F"s are arbitrary functions. This solution is a perturbative approximation to 
Vachaspatis' travelling waves 

For the solution with uj = and n = l,p = 0, the dispersion relation in the NLHE case is i? — fc^ = 0. The line of zero scalar 
field is deformed from the z-axis, w{t, z) = 0, to the helix w{t, z) — (,{k'^, k, 0) exp(— fc^i -I- ikz). Any initial disturbance of this form 
(for example a helical standing wave) will shrink down to the unperturbed straight-linear vortex. 

We have found some bound states for a; < 0. Let me explain the method first. The asymptotes at infinity are 



exp(-/MO , nexp(-|-/Mr) 

W_ — A r= h -D ■= , 

= C-M-VWYT2r] ^ ^eM+VWYT2r] 
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We want to find such w's that B and D vanish simultaneously. In Eqs.(^^ there is only one free parameter a. For each cu we can 
choose such a value of a that the coefficient D vanishes. Once D is removed the coefficient B remains in general nonzero. This 
procedure may be used to determine the function B{uj) = B{lu \ D = 0). We assume the function to be continuous. If there are 
two w's with opposite signs of B, there must be such a value of lo in between that B ^ 0. This is the way in which one can find a 
bound state solution to Eqs.(p^). 

For n — 1 and p = we have scanned a wide range of negative w's finding the function B to be always positive. More fruitful 
was the search for n = 2,p = 0. We have found a bound state for uj — —0.168 = —kg. The static modes which arise from this state 
are the same for all the considered models. The modes satisfy fc^ — fcg. The lines of zero scalar field are in general the roots of 
the complex polynomial w"^ — ^le**'"^ — ^26"**'°^. If only ^1 7^ 0, the solution is a static double helix, see the bottom part of Fig.l. 
Another characteristic solution is that with ^1 = ^2- This is a static double braid, see Fig. 2. The static double helices or braids 
exist only for the special characteristic wave-lenght Lq = 27r/fco- Perturbations with other wave-lenghts are no longer static. 

For the NLHE the dispersion relation is E = — kg. E is negative for fluctuations with wave-lenghts greater then Lq. The 
external radius of such double-helices or braids is growing with time like exp{^Et). Because of the string tension they must stabilise 
at some larger radius which we are not able to reach by our perturbative analysis. For short wave-lenght fluctuations E is positive 
and such double-helices or braids shrink down like exp{—^Et) to the unperturbed straight-linear vortex conflguration. 

The NLWE has a different dispersion relation, namely E"^ — k"^ = "^o- F'or generic values of positive k^ and E^ the double- helix 
at the bottom of Fig.l and the braid in Fig. 2 move up or down the z-axis with a phase velocity less then 1. Such modes should 
have been expected since they are just Lorenz boosted static helices or braids. More interesting are the modes with negative E^: 
—E^ = kg — k^. They exist for long wave-lenght fluctuations as compared to the critical wave-lenght Lq. Such double helices can 

decay or expand with a passage of time. There is also a combination for which the external radius evolves like \J sinh(V —E'^ t). 
This can describe a double-helix (braid) shrinking to the basic solution and then expanding but instantenously rotated around the 
z-axis by the right angle. 

In the limit of infinite wave-lenght we obtain a pair of parallel vortices. Vortices in the NLHE repel one another and the 
separation of their zeros grows like exp(ifcoO- If they started right from the unperturbed solution they would need infinite time to 
separate but any generic initial dipole fluctuation should substantially lower this time. For a similar pair in the NLWE, the most 
interesting time evolution is that in which the relative position evolves like •\/sinh(fcot). The positions of the zeros of the scalar field 
are the same as the zeros of the complex polynomial — sinh(fcof). For negative times vortices approach one another along the 
y-axis, at the time t=0 they concide and for later times they split but this time along the x-axis. This is the right angle scattering 
in the head-on collision of two vortices. This result can be generalised to symmetric collisions of n vortices to give the ^ scattering. 
A physical importance of this simple right-angle scattering has been realised by Manton 0. Its net effect for multivortex system's 
thermodynamics is that vortices behave as if they had finite cores - there is a net excluded area for a planar vortex liquid Q . 

3.3 Bound states for p = n 



We have seen, see Eq.(21), that the asymptotic equations for and u_ decouple. For p = n the two equations decouple for any 



r. Eqs.(|14D can be rewritten as 

LUU^ = — A„M+ + (3p — l)u+ , 

uju- = — A„u_ + (p — l)u_ . (25) 

What are the bound states of these stationary Schrodinger equations? As U- is concerned, we have not found any bound states for 
n = 1 and n = 2. The proof in Section 3.1 of Ref.|^ applies here. According to this argument there are no bound states for u_. 

The results for u+ are more interesting. For n=l there is one bound state with energy uji — 1.806. For n=2 the potential (3p— 1) 
is broader so that the bound state energy is lower and amounts to UJ2 — 1.613. It seems that the energy of the ground state is 
decreasing with increasing winding number. For sufficiently large n the next bound state is likely to appear. The existence of at 
least one bound state for any n can be proved following the lines in Section 4 of Ref . |^ . 

The two modes look qualitatively the same for n — I and n = 2. Once again let me begin with static modes. They satisfy 
k'^ — —LOn. The vortex exponentially broadens or tightens with z. An example is shown in the top part of Fig.l. Such solutions 
may seem at first sight to be simply diverging outside the perturbative regime and thus useless. However looking at Fig. 2 one can 
imagine that there is a surface of the superfluid somewhere above and our mode is the way in which the vortex mouth begins. 
Vortices prefer to be wider at the surface then in the bulk. 

Now about the time-dependent solutions. In the NLHE the dispersion relation is E = Un + k'^ ■ One can perturb the vortex 
changing its core radius but the distortion will shrink down to the unperturbed solution like exp—Et. It shrinks the faster the 
shorter is the wave-lenght of the fluctuation. In the NLWE, where the dispersion relation is E^ ~ k^ — ojn, such a z- independent 
distortion oscillates in time. Something like a "breather" or "pulson" state forms. Making the last solution z-dependent, with k^ > 0, 
we obtain waves travelling along the vortex with a phase velocity greater then 1. 

There are also solutions of Eqs.(|25|) just on the verge between bound states and scattering states. One of them is u_ — f{r) for 
uj = 0, where /(r) is the moduli of the background solution. Another solution is u+ = f{\/3r) for oj — 2. These solutions would exist 
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for any background, they have nothing to do with the vortex. The first solution is a zero mode due to a global gauge transformation 
4> — > tpe^^° . The second one would be u+ = 1 for a uniform background p = 1. It is just an uniform massive oscillation of the scalar 
field's modulus around its equilibrium value equal to 1. 

3.4 Bound states for p < 

As discussed in Section 3 such modes describe splitting of the n-vortex into n unit vortices and in addition | p \ vortex-antivortex 
pairs. Such a decay is in principle possible because the energy of n-vortex tt?) is higher then the energy of n unit vortices (~ n). 
The extra energy can be used to create vortex-antivortex pairs. The maximal number of such pairs has been estimated to be at 
best I p 1= n — 1, see Q. 

For n=l and n=2 the result is negative. Double- vortex can be split into 2 unit vortices but there is not enough energy to create 
a. V — V pair. Such a decay turns out to be possible for n = 3,j3 = —1. There is one antivortex left at the origin and 4 uniformly 
split vortices. The corresponding eigenvalue is a; = —0.103. For the NLHE the mode is a diffusive splitting with creation of one 
V — V pair. This way of splitting has been observed in a direct numerical simulation, see the Figure 3 in |6j. z-dependence can 
be introduced giving rise to vortex 4-helix around a straight-linear antivortex or 4-braid with nodal points on the antivortex. For 
fc^ = —Lj the helices and braids are static and for larger fc^ they are damped. For the NLWE there are 4-helices and 4-braids at the 
critical wavelength. At other wavelengths they are travelling waves. For parallel vortices the mode describes head-on collision of 3 
vortices on an antivortex center resulting in ^ scattering. 

4 Fluctuations in the nonlinear Schrodinger equation regime 

The discussion of the NLSE has been postponed until this section. The NLSE can be obtained from the NLHE by a formal Wick 
rotation t = ir 

i^riJ^-V^^p + i^p*^p-l)^|J . (26) 

However it does not mean we can do the same with solutions of these equations. For example the eigenvalue problem for the 
fluctuations around a vortex background in the NLSE 

^{E, k,p)fo[e'^*+'''''e'P%i{r) + e-'-^*-*'=^e^(2"-p)%2(r)] (27) 

takes a slightly different form 

{-E + + 2p-l- Ap)ui + pu2 = , 

i + E+k'^ +2p-l- A2n-p)u2+ pui^O (28) 

than in the NLHE. Note that the equations differ by a sign in front of E. This makes the eigenvalue problem different from that 
for the NLHE and NLWE, compare with Eqs.(|l3|). This is where the analytical continuation in the complex time plane fails on the 
perturbative level. 

Let us take a closer look at asymptotic properties of the solutions to the eigenvalue problem (^) . The asymptotes at the origin 
are the same as in the NLHE case, see Eq.(p^. As E and fc^ can not be combined now into just one variable iu — E — k'^, the 
asymptotic behavior at infinity will depend on both E and fc^ independently. In this limit the Eqs.(|2^) become 

Aiti + {E -k'^ ~ l)ui -U2 = , 

Au2 - -iti + {-E -k^ - l)u2 = . (29) 
These equations can be diagonalized. The eigenvalues turn out to be 



-(l + fc2) - v/l + S2 . (30) 

For a positive eigenvalue we obtain a Bessel function-like asymptote which does not converge fast enough to be acceptable. The 
two eigenvalues are negative for (i?, fc^) belonging to the area defined hy k^ > — \ + \/l + E^. In this "convergence area" there are 
two exponentially decaying asymptotes and two exponentially divergent ones. For any {E, fc^) one of the divergent asymptotes can 
be removed with an appropriate choice of the constant a in Eq. (p^) . The other one may happen to vanish along some lines in the 
{E, fc2)-plane. These distinguished lines are just the dispersion lines we are looking for. 

A straightforward but laborious way of solving the problem ( p8| ) would be to scan the {E, fc2)-plane in search of bound state 
solutions. Let us check first what could we learn from analytical continuation of the NLHE's bound states to the NLSE. 

Let us consider the bound state of the NLHE 
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where the functions ui(r),U2(r) satisfy Eqs.([l3|). The Wick rotated configuration takes the form 

c~/ —lEtr ikz ip9 /\ . —ikz i(2n — p)6 /_M /"ooA 

dip = e [e ui(r) + e e ^ ^' U2(r)\. [32) 

This configuration is a solution of the hnearized NLSE but only at t = 0. Later on the fluctuation 6ijj must deviate from the 
continued mode Sij) to 6ip = Sip + cj). We look for the deviation in the form 

The profile functions w's and w's must satisfy the following sets of inhomogenous differential equations 

{+E + P + 2p-l- Ap)wi + pw2 = -2pu2 , 

i-E + + 2p-\- A2n-p)w2 + pwi = 2pui (34) 

and 

(-E + fc2 + 2p - 1 - Ap)ui + pv2 = 2pu2 , 

i+E + fc^ + 2p - 1 - A2n-p)v2 + pvi = -2pui . (35) 

The R.H.S.'s of the above inhomogenous equations are regular sources which vanish exponentially at infinity. Close to r = the 
w's and w's look like 

(36) 



■wi{r) f 




W2{r) f 




vi{r) R 




V2{r) R 





with w'^'s and w'^'s being constants which have to be choosen so as to make the solutions convergent at infinity. This time the 
equations are inhomogenous so this time the constants can not be rescaled to remove a half of them. For each set of equations 
"35]) there are two adjustable parameters, (wj'jW^) and (i^PjU^) respectively. As to the asymptotic properties at infinity, it is 



enough to note that the differential operators in both sets (34 35|) are the same as in the Eqs.(|28|). If the eigenvalues (i?, k"^) belong 
to the convergence area fc^ > — 1 + vT+~E^ than, for each set of equations ( |3^ , |35| ), there are two divergent asymptotes at infinity. 
These two asymptotes can be removed with the two free parameters in the asymptotes close to the origin (^6|). 

Thus it turns out that is a quickly convergent function provided that k"^ > — 1 + y/l + E^. For (E,k'^) belonging to this 
convergence area the NLHE mode 6ip when analytically continued to the NLSE gives rise to the mode 6ip — 6ip + cj). This means 
that the parts of dispersion lines which belong to the convergence area can be analytically continued from the NLHE to the NLSE. 
We have confirmed this observation looking directly for the solutions of Eqs.(p8|). Let us then discuss the most characteristic 
solutions of the NLSE. 

For n=2 and p=0 the NLHE's dispersion relation is E — k^ — k^. The dispersion relation remains unchanged for the NLSE 
modes obtained by mapping NLHE's solutions. For short wave-lenght double helices and braids (Fig. 1, Fig. 2) the energy is positive. 
They can move up or down the z axis as travelling waves. For fc^ = kg they become static. In the range of k^ from fcg to 
i(l - v^l - 4A:g) there are once again travelling waves. Now the problem arises if we can continue the dispersion line E = k^ — k^ 
ouside the convergence area, in particular to the point {E — —kg, fc^ = 0) which corresponds to a pair of parallel vortices rotating 
anticlockwise around their common center of mass with angular velocity ^k^. Outside the convergence area the long-range Bessel 
function-like asymptotes are unavoidable. A well defined NLHE's mode when Wick rotated to the NLSE developes a long range 
deviation 4>. Thus it seems that a rotating parallel pair of vortices dissipates energy radiating sound waves. As the potential between 
vortices is repulsive their mutual distance should be growing with time. 

For the excitations (n = l,p = 0) of a single vortex the dispersion relation is -E = fc^. This line remains in the convergence area 
for all k^ > 0. The travelling waves on the single vortex move up or down the z-axis. The phase velocity is falling down to zero as 
the wave-lenght tends to infinity. These excitations can be identified with the experimentally observed Kelvin modes 

The dispersion relation for the breather states {p = n) is E = Un + fc^ with ujn positive but smaller then 2. For n = 1, 2 the 
whole dispersion line lies outside the convergence area. It may happen that for some large enough n, w„ < 1. Then there would 
exist a finite part of the dispersion line belonging to the convergence area. In any other case the breather modes must radiate. For 
A:^ > they travel up or down the vortex with some phase velocity dependent on the wave-lenght. In the z-indcpendent case fc^ = 
there is a breather state like in the NLWE - vortex width is oscillating in time around the equilibrium value. For fc^ < there exists 
an exceptional point {E = 0, fc^ = a;„) where we have a static vortex throat common to all the considered equations. 

In the n = 3,p — —1 case 4 vortex lines uniformly split from the origin. For k^ large enough, the 4-helix or the 4-braid wave is 
travelling along the central antivortex axis. The point on the dispersion line with fc^ = lies outside the convergence area thus the 
mode with parallel vortices rotating around the central antivortex is at best a radiative one. 

The fluctuations are similar to those in the NLWE. The difference is that the dispersion relations for the modes in the NLSE 
are a nonrelativistic limit of the dispersion relations for corresponding modes in the NLWE. Some of the NLSE modes, if they exist 
at all, are unstable against decay by radiation of sound waves. 
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5 Mixed Schrodinger-difFusive dynamics 

So far we have considered only the NLHE and NLSE. In general we can replace t in the NLHE ( p^ by tex.p{ij) to obtain the 
mixed Schrodinger-diffusion equation 

— i/i^t COS7 + zi/'.t sin7 = — A-)/' — (1— I ?/> P)?/" ■ (37) 
7 can vary from to ^ . The fluctuation modes now take the form 

fc,p)/oe-^*'^°'^^[e-*^*™'^e''^'^e'P%i(r) + e+'-^*™''e-''=^e'(2"-P)%2(r)], (38) 

where this time ui, U2 are complex. To answer the question if the exponentially localised NLHE's modes can be mapped into modes 
of the intermediate equation we have first to find its convergence area. 
The homogenous equations which ui, U2 have to satisfy are 

[Ap - fc2 + (1 _ 2p) + E]ui -pu2=Q , 

[A2„-p -k^ + 2p) + £;exp(-i7)]u2 - pui = . (39) 
The real part of the eigenvalues of the matrix operator (for very large r) is 



- [l + k^) + Ecos^ + Re^l- E'^siT? -ie-'^''^ . (40) 
The convergence area is a set of all such {E, k^) that 



k"^ > -1 + Ecos^-/+\ Re^Jl- E^sm^-/e-^'f \ , (41) 

where the real parts of both eigenvalues are negative. In this area, by similar arguments as in the previous section, a deviation (p 
from a continued NLHE's solution can be always made convergent. 

The dispersion relations are once again formally the same as for the NLHE or the NLSE but as for the NLHE the sign of E 
determines whether a given mode is growing or decaying with time. 

A pair of parallel vortices not only rotates around the common center of mass with the angular velocity ifcp sin7 but also 
their mutual distance is growing with time. The zeros of the scalar field move along spiral lines. There is such a 57 > that for 
7 < ^ — (57 the dispersion relation E — k'^ — kg for the {n — 2,p — Q) mode can be continued from large positive k^ down to the 
axis fc^ = and a little bit below. Thus a pair of parallel vortices rotating one around another is radiating sound waves in the 
NLSE regime but a small dissipation makes them an exponentially localised solution. The external radius of the long wave-lenght 
double-helix (Fig.l) or the braid (Fig. 2) is also growing with time until it stabilises thanks to the string tension. If its wave-lenght 
is short, it prefers to shrink down to the coincident two- vortex configuration. 

For single vortex modes {n — l,p = 0) the dispersion relation E = k^ remains in the convergence area for all positive k^. A 
single-helix gradually shrinks down to the straight linear vortex. Fluctuations in the vortex width, the same as in the NLSE case, 
remain outside the convergence area. 



6 Summary 

We have obtained a wide class of perturbative solutions like helices, double-helices or braids. The question arises whether such 
configurations really correspond to some exact solutions. We suppose the answer to be yes. Configurations like helices (Kelvin 
waves) or double helices were analysed within the models for widely separated vortices (Kelvin waves are exact solutions to the 
NLWE [^). The braids can not be described by such models because even if their amplitude is large vortices have to cross one 
another at the nodal points. In the Bogomol'nyi limit of the Abelian Higgs model some extra symmetries of the model enable an 
exact construction of double-helices and braids 

The analysis was substantially simplified by various symmetries which connect some of the considered equations. The perturba- 
tive calculations in the NLWE and the NLHE regimes lead to the same stationary eigenvalue problem. The perturbative solutions 
to the whole family of nonrelativistic equations interpolating between the NLHE and the NLSE were constructed from the NLHE's 
solutions. This construction shows that although the analytical continuation of solutions in the complex time plain fails one can 
still continue the dispersion relation, provided that dispersion line belongs to the convergence area. The calculations in the NLHE 
were substantially simpler then they would be in the corresponding problem for the NLSE. The lesson is that it is convenient to 
map such problems to the NLHE and later on continue the obtained NLHE's dispersion lines to the equation under consideration. 

Note added. A month after the first version of this paper the preprint by Goodband and Hindmarsh | pT[ appeared. The papers 
overlap in the analysis on the NLWE. 
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7 Figure captions 

Fig.l. Superposition of two modes on a double vortex. The bottom part is the double helix. In the top part the vortex mouth opens. 
On the plotted surface the scalar field's modulus is equal to i. 

Fig. 2. The double braid. The double vortex splits into two unit vortex cosine waves of the same polarisation. At each nodal 
point the polarisation plane turns by the right angle. For n-braids the polarisation would turn by the angle - at each nodal point. 
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